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The mutual drag in strongly interacting two-component superfluids in optical lattices is dis- 
cussed. Two competing drag mechanisms are the vacancy-assisted motion and proximity to the 
quasi-molecular state, in which an integer number q of atoms (or holes) of one component might be 
bound to one atom (or hole) of the other component. Then the lowest energy topological excitation 
(vortex or persistent current) becomes a composite object consisting of q circulation quanta of one 
component and one circulation of the other. In the SQUID-type geometry, the value of q can become 
fractional. These topological complexes can be detected by absorptive imaging. We present both 
the mean field and Monte Carlo results. The drag effects in optical lattices are drastically different 
from the Galilean invariant Andreev-Bashkin effect in liquid helium. 



Multi-component quantum mixtures in optical lattice 
(OL) are a source of new and rich many-body physics. 
The one-component bosons in OL have been exhaus- 
tively studied theoretically and experimentally 0, 
especially in the context of the quantum phase transi- 
tion between superfluid (SF) and Mott insulator. The 
study of multi-component boson systems in OL has just 
began. Theoretical investigations predict variety of new 
quantum phases with unusual properties 0,0,IEIES@- 
Two interesting recent examples include topological ex- 
citations - vortices and persistent currents with non- 
standard winding properties in two-component superflu- 
ids (2SF) HQ. 

Crucial, but largely unaddressed effect is the impact 
of strong interaction on properties of superfluid phases 
where each component -0a has its finite expectation value 
(ipa)- Interesting manifestation of the strong interac- 
tion is the inter-component drag similar to the Andreev- 
Bashkin effect in 2SF helium mixtures. In general, 
the drag effect between non-convertible species at zero 
temperature is represented by the cross-terms in the ex- 
pansion of the ground state energy in terms of small gra- 
dients of the superfluid phases V<^ a , a = 1, 2, 



SE = 
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with p ao standing for the superfluid stiffnesses. The 
cross-term p\ 2 is responsible for the drag. It is due to 
interaction effects and is not confined to some particu- 
lar term in the full microscopic many-body Hamiltonian. 
Depending on its sign, this term describes either a mutual 
unidirectional flow or a counterflow of the components. 
The Galilean invariance argument, often attributed to 
Landau, imposes two constraints on p ao . These con- 
straints are responsible for the Andreev-Bashkin effect in 
superfluid mixtures of liquid helium isotopes in which p 12 
is uniquely related to the ratio of bare TOi and effective 
to* masses of minority atoms in the host superfluid of the 
majority component. The Galilean transformation to a 
frame moving with velocity V requires that the phase of 
each component changes as tpi — > (pi—(mi/H) V ■ r, where 



mi are the bare masses. The energy density Q trans- 
forms as SE -> 6E — P V, where P/H= N x \Iodi+N 2 Vld 2 
is the momentum density expressed in terms of particle 
densities Aq^ of each component. This yields 

pumx + p\ 2 m 2 = Ni, pi2TOi + p22m 2 = N 2 . (2) 

Introducing effective masses m 1; m 2 as p\\ = 
Ni/ ml, p 22 = A 2 /m*., we reproduce the result 
P12 = {Ni/m 2 )(l — mi /to*) as well as the relation 
Ni(mi — ml)mi/ml — N 2 (m 2 — m^rn^/rn^- In other 
words, conservation of the total momentum requires that 
the difference of the bare and effective masses is compen- 
sated by the flux of the other component. Note that 
P12 > since to* 2 > m\ y2 . 

In the case of strong mass renormalization, (to* / m i ) — 
1 > 1, quite spectacular effects should be expected [lCj 
from the topological excitations - vortices. Specifically, 
the lowest energy single-circulation vortex of the majority 
component (p 22 3> pn) should carry several circulation 
quanta q = 1,2,.. of minority component. The equi- 
librium value of q is obtained by minimizing the factor 
m 2 q 2 + 2 (ml — toi) q in the energy of the vortex com- 
plex (or persistent current). These q+1 vortex complexes 
exhibit transformations with respect to the value of q de- 
pending on external conditions that determine the value 
of to*. If the interaction is weak, pi 2 can be calculated 
as an expansion in the gas parameter pd| . 

In this paper, we address the drag effect in a lattice 2SF 
in strongly interacting limit, and show that it is radically 
different from the Galilean-invariant case. The lattice 
plays a central role in violating the relation between 
pi 2 and TOi/to* (and the constraints (|2J)). We also argue 
that the value of q is affected by proximity of the 2SF to 
the quasi-molecular phase. 

In OL, in contrast to the Galilean-invariant system, the 
lattice provides a preferred reference frame, so that the 
(hydrodynamic) properties of the two-component mix- 
ture are determined not by the relative velocity of com- 
ponents but by their individual velocities with respect 
to the lattice. Furthermore, the effective mass in OL is 
formed largely by the width and depth of laser-generated 
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potential wells rather than by a trailing cloud of the sec- 
ond component. Another crucial difference is that in 
OL number of vacancies is a conserved quantity. Be- 
low we perform the mean field and Monte Carlo analy- 
sis of the mutual drag in 2SF in three different physical 
situations: a soft-core system close to molecular conden- 
sation, a hard core system with finite intercomponent 
exchanges, and a hard core system with vacancy-assisted 
motion without the intercomponent exchanges. 

Drag due to proximity to the quasi-molecular state. 
Here we discuss a generic mechanism leading to the 
q + 1-topological complexes in the 2SF. Strong drag ef- 
fect occurs if a two-component boson system is close 
to a transition into the quasi-molecular state in which 
the only broken symmetry has the order parameter 

$ g - exp(^ (<j) ) ~ (iptipz) + ( or (V-i^l 9 ) ^ 0). In 
pure molecular state with undefined individual phases 
ipi t 2 (that is, ("01,2) = 0), the phase-gradient energy is 
given by the molecular superfluid phase tp( q > as SE = 
f dxp q (Vtp t - q ' > ) 2 /2, with p q being molecular superfluid 
stiffness. The molecular order parameter persists in the 
2SF phase so that the additional broken U(l) symmetry 
emerges continuously [||. The two phases y?x,2 become 
well defined in the 2SF state with the molecular phase 
being locked as 

ip (q *> = (fix + qip 2 . (3) 

This locking can be understood as a consequence of vir- 
tual processes of transformation of a (g+l)-molecule into 
q B-atoms and one A-atom. The corresponding contri- 
bution to the energy functional is AE ~ j dx^ g ip*ip2 9 + 
H.c. This term (c/. the diatomic molecules with q = 1 
[13) ensures the relation (JSJ in the longwave limit. Then 
the energy Q becomes 

6E = J dx[^(V(^+^ 2 )) 2 (4) 

+ y(V^) 2 + ^(V^ 2 ) 2 + pi 2 V^V^ 2 ], 

with p'lj continuously changing from zero in the molec- 
ular phase to some finite values in the 2SF phase. It is 
important that the molecular stiffness p q is not a crit- 
ical property of the system - it does not change while 
crossing the phase boundary. Thus, at least close to the 
phase boundary, minimization of the vortex energy gives 
ipi = —q(f2, that is, the q + I vortex. In reality, the re- 
lations \p' ab \ <C |/t>x2 1 ~ Pi, 2 can hold quite far from the 
phase boundary. This implies that the q + 1 topological 
excitation exist deep in the 2SF phase. We demonstrate 
this numerically for q = 1 (see Fig.l below). 

It is convenient to introduce the drag coefficient k as 
a ratio k = pvxj p\\ of the cross-stiffness to the smallest 
diagonal stiffness, p\\ < p2 2 . Then, as the minimiza- 
tion of the energy JIJ shows, when |fc| > 0.5, a vortex of 
the dominant component can lower its energy if it car- 
ries the circulation of the other component q = ±1. In 
symmetric case (pxi = P2 2 ), the integer q closest to k 



determines the q+1 vortex (or persistent current) as the 
minimal topological excitation. It is important to note 
that even small |fc| causes attraction between either vor- 
tices of equal circulations (k < 0) or between vortex and 
anti-vortex (k > 0) in different components, so that if 
both exist they will form a complex. Crossing the bound- 
ary | Ac | > 0.5 has strong impact on mechanisms of vortex 
creation and stability. For example, stirring the compo- 
nent with the largest stiffness (P22) above the threshold 
will cause creation of the complex instead of a single vor- 
tex of the stirred component. Also, a single vortex of the 
component 2 becomes unstable with respect to inducing 
creation of vortex of the other component. 

The Hubbard lattice model with molecular phases, 

H= ^2 h t a a a,i a aj + H.c] + [U a , a <n ai n a > ,j] , 

(5) 

has been extensively studied analytically 0, @ and nu- 
merically 0, 0. Here U ata i is the interactions matrix, 
t a describes the nearest-neighbor jumps of component 
a; a ai , a a ,j are the construction bosonic operators, and 
n>ai = a a i a a i are the on-site occupancies. As dis- 
cussed in Refs.j^Q, the quasi-molecular phase (U12 < 0), 
namely, the paired superfluid, is in many respects isomor- 
phic to the super-counterfluid state ({7x2 > 0) |4| . Both 
states can undergo second order phase transition into the 
2SF phase so that the order parameter ^^x = (tpi ^2) 
(or $ g= x = (V'xV'J)) remains finite and robust. Obvi- 
ously, in the 2SF phase, the q = 1 composite vortices 
are the lowest topological excitations. As pointed out in 
ref.[l3j, the Hamiltonian (JSJ also allows molecular phases 
with arbitrary integer value of q. This issue, though, re- 
quires separate analysis. 

Hard core limit U a b — > 00 of the Hamiltonian (JSJ) . This 
limit can exhibit quite interesting physics of strong quan- 
tum fluctuations even far from any phase transition [j|. 
Obviously, when Na + Nb = 1 (Na, Nb are the aver- 
age on-site occupancies of the species A, B), the system 
in the hard core limit (HC) is the Mott insulator. Its 
ground state is degenerate with respect to possible per- 
mutations of bosons A and B. This degeneracy, which is 
a consequence of the HC approximation, is lifted by any 
infinitely small inter-component exchanges. Accordingly, 
the two-component HC model should be considered as a 
limit of the model in which the inter-component inter- 
action Vi n t — {7x2 is finite and increasing. In contrast 
to free space, increasing Vi n t leads to decrease of all su- 
perfluid stiffnesses because all transport is suppressed as 
~ t\t2/Vint. Furthermore, in the limit Vmt — * co all stiff- 
nesses are equal in magnitude. This is clearly at variance 
with the free space constraints (J2J which prohibit uniform 
decrease of all stiffnesses at fixed densities. 

The two-component Hamiltonian with residual soft- 
core inter-component repulsion is represented in terms 
of the HC construction operators aj, (H and b], bi with 
Pauli commutation relations for the A and B components 
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H = Hl4 a j ~ l A hi + Hx -] + J2 V int4 a i b i b i ( 6 ) 
<ij> i 

with summation < ij > over the nearest-neighbor sites. 
At total filling 1, this Hamiltonian has two phases - 2SF, 
where both SF order parameters are defined, and super- 
counterfluid (SCF), where the only SF order is observed 
in (ciib'j). Transition between these two phases is con- 
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FIG. 1: The drag coefficient k for the J-current analog of the 
Hamiltonian © as a function of the relative interaction, with 
the value Vi n t/V c — 1 corresponding to the 2SF-SCF phase 
transition. The horizontal dashed lines indicate a domain 
where 1 + q- vortex complex with q = 1 has lower energy than 
any single circulation vortex. Error bars are much smaller 
then symbol sizes. Solid line is the eye guide. 

tinuous in the universality class U(l) |(| and occurs in 
the symmetric case t\ = t% = t at some value Vint = V c , 
V c /t ~ 1. As discussed above, the drag effect is strong in 
the 2SF phase even far from the transition. We proved 
this by performing the Worm algorithm Monte Carlo 
simulations of the two-color J-current model @, 0, El at 
zero temperature on a 2D square lattice. This model is 
a discrete-time grand-canonical analog of the Hamilto- 
nian 10 with the hard-core constraints. The stiffnesses 
were determined from the statistics of the winding num- 
bers similarly to Refs. ((JO, El- The SCF phase was 
identified by observing pu = P22 = —pu- The negative 
value of pi2 is due to counterflow of the components - 
each winding of A-worldline is accompanied by opposite 
winding of B-worldline. In Fig.l, the drag coefficient k is 
plotted as a function of the relative interaction strength. 
As can be seen, the domain l/2<|fc|<lin the 2SF (be- 
tween the dashed lines), where the composite 1 +q- vortex 
with q — 1 has lower energy than any single vortex, is 
not restricted to the vicinity of the critical point V c but 
occupies about half of the phase diagram. Here p\2 < 0, 
indicating that both components participate in the coun- 
terflow even in the 2SF state. 



Vacancy assisted drag. If the total filling is different 
from 1, system is always in 2SF phase at T = 0. In this 
case, another mechanism contributes to the drag - the 
vacancy assisted transport. Atoms tunnel to the unoccu- 
pied sites (vacancies) much faster than the rate of the 
A-B exchange with large Vi„t- The vacancies stimulate 
mass flow in one direction and move in the opposite one. 
As a result, both components A and B move in one direc- 
tion, which means that pi2 > 0. This situation implies 
crossover when pi2 changes sign at some special point 
with no drag, p\2 = 0. Since no symmetry change takes 
place, this is not a phase transition. The crossover from 
k < to k > takes place as Vi n t increases at fixed 
number of vacancies. 

Note that the drag coefficient k must increase when the 
number of vacancies x v = 1 — Na — Nb decreases. This 
counterintuitive result stems from the nature of vacan- 
cies. In one component case, conservation of the number 
of vacancies Ny makes them similar to particles. The 
HC limit links the flow of vacancies with the opposite 
flow of atoms. In the two-component case, the situa- 
tion is similar with one crucial difference - a vacancy is 
not uniquely associated with a particular sort of atoms. 
Thus, motion of a single vacancy through a lattice in one 
direction leads to flows of both components in the oppo- 
site direction. This implies strong drag with positive k. 
When x v increases, system becomes more like a low den- 
sity and, thus, weakly interacting mixture of two sorts of 
atoms with correspondingly small k. 

To analyze the mutual drag and the possibility of com- 
plex vortices in the vacancy dominated regime we modi- 
fied the HC model by imposing the additional constraint 
aibi = on @ and introducing the chemical potentials 
term —p,\NA — P2^b for each component to keep con- 
trol of the filling factors. As discussed in Ref.Q, this 
limit can exhibit long range phase separation as well as 
short scale fluctuative phase separation corresponding to 
minority particles acquiring large cloud of vacancies. 

If 1 — Nb > Na, it is convenient to introduce a de- 
scription in which the vacuum corresponds to all sites 
filled by B particles. Then, the number n = I — Nb 
of B holes is shared between Na atoms and remaining 
x v = n — Na > vacancies. In the limit Na *C Nb ~ 1 
transport of vacancies can be considered as transport of 
B holes with the effective Hamiltonian 

<ij> 

where v\, Vi are the Pauli operators for B holes. In order 
to describe the mutual drag within the mean field ap- 
proximation, one should replace the field operators a, v 
by the functions a = y/jTfexp (*</?i), V — v / x^'exp (—i(p2) 
with the slowly varying phases and perform the gradient 
expansion. [The minus in front of iip2 indicates that flow 
of holes and actual flow of mass are opposite] . This au- 
tomatically generates the term ~ tiXix v (V(ipi + 1P2)) 2 
in effective energy from the first term in Hamiltonian 
0. Obviously, the ratio of the stiffnesses becomes 
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k = Pw/pu = 1 which corresponds to positive cross- 
term typical for the vacancy assisted transport meaning 
that the mean field captures well the physics of the va- 
cancy assisted transport. However, the prediction k = 1 
and, therefore, q = 1, is not supported numerically. 

We have performed Worm algorithm 0] Monte Carlo 
simulations of the two-color J-current model [(J 0, El at 
zero temperature in 2D square lattice in the HC limit 
with partial filling. This model is similar to described 
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FIG. 2: The drag coefficient k for the J-model analog of 
the Hamiltonian JHJ in the limit Vint -* oo as a function 
of concentration of vacancies x v for symmetric case: Na = 
Nb, t a = tt- Error bars are shown for all points. Solid line is 
the eye guide. 

above 10 with an additional requirement of no double 
occupancy. We have found that (Fig. 2), for x v > 0.15, 
k < 0.42 ± 0.02, and, thus, no topological complexes can 
exist as the lowest energy topological excitations in this 
regime. At this point we do not have a simple explana- 
tion for this variance between numerical and mean field 



results. Most likely, the mean field result is not appli- 
cable for large x v in the symmetric mixture when the 
vacancies cannot be uniquely identified with the holes in 
the majority component. 

As the number of vacancies is tuned to become x v < 
0.15, all stiffnesses exhibit large error bars which can be 
attributed to the regime of strong quantum fluctuations 
[B| associated with the degeneracy of the ground state in 
the HC limit. The precise nature of this effect requires 
separate analysis. For finite Vint-, depending on Na, Nb, 
the ground state can exhibit various types of ordering in- 
cluding the checkerboard insulator Q. Then, decreasing 
x v at Na = Nb — ► 0.5 will result in the first order phase 
transition with strong fluctuations, similar to those in 
Fig. 2, due to the domain formation. 

Fractional q. In the case of finite drag with |fc| < 0.5 
fractional phase circulation q — k can be observed when 
persistent current is interrupted by a Josephson junction 
which lifts the requirement of the integer of 2tt windings 
by creating the phase jump across the junction. Then, 
phase winding is determined solely by the minimization 
of energy. 

Detection. Finally, the (q + 1) -vortex complexes can 
be observed by absorptive imaging technique similar to 
imaging of vortices in one-component Bose-Einstein con- 
densates |17| . Typical pattern should include extra q 
fringes in one component. 

In summary, we explored generic mechanisms of drag 
effect in quantum bosonic mixtures in optical lattice with 
hard and soft core interaction. Strong mutual drag can 
result in composite topological structures. The drag in 
lattice is not controlled by particle effective masses. The 
simplest mean field approximation does not adequately 
describe the strong drag. 
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